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Lenses that can collect the perfect image of an object must restore propagative 
and evanescent waves. However, for efficient information transfer, e.g., in com¬ 
pressed sensing, it is often desirable to detect only the fast spatial variations 
of the wave field (carried by evanescent waves), as the one created by edges 
or small details. Image processing edge detection algorithms perform such 
operation but they add time and complexity to the imaging process. Here, we 
present a new subwavelength approach that generates an image of only those 
components of the acoustic held that are equal to or smaller than the operating 
wavelength. The proposed technique converts evanescent waves into propaga¬ 
tive waves exciting trapped resonances in a waveguide, and it uses periodicity 
to attenuate the propagative components. This approach achieves resolutions 
about an order of magnitude smaller than the operating wavelength and makes 
it possible to visualize independently edges aligned along different directions. 
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Edge detection is an essential numerical tool in image processing that finds application 
in several areas of science and technology. In medical imaging DEI. non-destructive 
testing M, and computer vision [5], edge detection plays an important role, since it 
enables extracting the meaningful information from an image and it reduces the amount 
of data to be processed. The basic idea behind this technique is to high-pass filter the 
image to remove the low spatial frequencies. Close to the edges of an object illuminated by 
a monochromatic wave, the wave field is dominated by evanescent waves, i.e., waves with 
spatial oscillations faster than the operating wavelength. A lens capable of generating an 
image using only evanescent waves would visualize the edges or small details of an object, 
essentially extracting only the key information contained in the image. 

There exist different ways to detect evanescent waves, which overcome the classical 
diffraction limit [6] of conventional imaging devices. Approaches based on superlenses 0 - 
GDj) and hyperlenses (HUSH! or time reversal techniques [221 fZo] allow restoring evanescent 
waves, providing a detailed picture of the imaged scene. However, since propagative waves 
that carry low spatial frequencies are also used to form the image, such techniques cannot 
be used to visualise only fast wave field variations. 

Here, we present a new imaging technique for acoustic waves, based on trapped res¬ 
onances in irregular waveguides, that provides an image of only spatial variations of the 
acoustic field that are equal to, or smaller than, the operating wavelength. This approach 
provides sharp images of the edge of an object, with resolution up to A/7.6 (A is the 
operating wavelength). Moreover, the technique allows visualizing edges aligned along a 
given direction independently. These features arise from the fact that the plane mode ex¬ 
citation is not involved in the restoration of the evanescent components. This is different 
from other imaging approaches that use acoustic resonances, which rely on transporting 
propagative and evanescent information mm- 
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We design waveguides with periodic, symmetric modulations of the cross-section, as 
shown in Figures 1A and IB. The waveguides have two different square cross-sections; a 
narrow section, s, and a wide section, S, with dimensions s = w x w, and S = W x IV. 
The axial lengths of the narrow and wide segments are denoted respectively by l and L. 
The symmetric variations of the cross-section generate trapped resonances (TRs) in the 
waveguide [26], [27], which appear at frequencies slightly above the cutoff frequency of the 
first antisymmetric waveguide mode. More details about the origin of these resonances are 
given in the Supplementary Section SI. These resonances are antisymmetric with respect 
to the longitudinal axis (see Supplementary Figure SI), meaning that they can only be 
excited by antisymmetric waveguide modes (some of those modes are shown in Figure 
1C). The excitation of the TRs induces a strong coupling between higher-order modes, 
including the evanescent ones, making it possible to tunnel subwavelength information 
through the device. The periodicity induces a bandgap for the plane mode, which avoids 
the transmission of components with small perpendicular wavenumber. By choosing the 
geometrical parameters of the waveguide, it is possible to make the plane mode band gap 
and the TRs coincide in frequency, creating a spectral band in which only waves with 
large perpendicular wavenumber are transmitted. In the particular case considered here, 
the geometrical parameters are w = 7.5 mm, W = 2.25 mm, l = 3 mm and L — 15 mm. 
Figure ID illustrates the principle of operation of the waveguie: low spatial frequencies, 
carried by propagative waves (blue sinusoidal lines), are converted into evanescent waves 
(red decaying lines). High spatial frequencies, carried by evanescent waves, are converted 
into propagative waves. As a result, a picture of only the edges of the imaged object is 
created. 

The acoustic field inside the waveguide is described by the 3D wave equation 
(V 2 + k 2 )p(x,y, z) = 0, where k = lo/cq is the wavenumber, ui is the angular frequency 
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Figure 1: Imaging device. (A) Schematic of the device, consisting of a waveguide with 
periodic modulations of the cross-section. The frontal waveguide wall is not shown in 
order to expose the internal structure. The holes pointed by green arrows indicate the 
microphones position in the experimental study. (B) Experimental realization using 3D 
printing. (C) Eigenfunctions of modes (0,0), (1,0), (3,0), (5,0) and (7,0) and their 
position in the frequency axis. At the trapped resonance frequencies (black dot) only 
modes (0, 0) and (1, 0) are propagative. (D) Basic illustration of the operating principle 
of the waveguide. 

and c 0 is the speed of sound in air, taken here as c 0 = 343 m/s (the time dependence 
exp(— jut) is omitted). Assuming rigid boundary conditions, the solution p(x,y,z ) can 
be expanded on the modal basis of the waveguide cross-section as 

OO 

P(x,y,z)= ^ A (m,n)(x)4>(m,n)(y, Z), ( 1 ) 

m,n=0 

where (?/, z) are the eigenfunctions, A^n^x) is the modal amplitude as a func¬ 

tion of x , and the couple (m, n) indicates the number of vertical (m) and horizontal (n) 
nodal lines. The propagation of each mode is determined by its longitudinal wavenumber, 
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f3(m,n) — with = [(mn/W) 2 + {u'K/W) 2 } 1//2 the transverse wavenum¬ 
ber. For a given frequency, these modes are propagative if k > or evanescent if 

k )• 

Using the mode-matching technique (see Supplementary Section S2), we have cal¬ 
culated the transmission matrix T, A 1 = T A 1 , where vectors A 1 and A T contain the 
incident and transmitted modal amplitudes. In Figure 2A we show the plane mode 
transmission coefficient, T(o,o),(o,o), and the term corresponding to the first antisymmetric 
mode, T , ( li o)(i,o), in the frequency band [0,12.5] kHz. The curve corresponding to mode 
(0, 0) shows a wide band gap between approximately 3.7 kHz and 11 kHz. In the absence 
of TRs, the transmission coefficients of any higher order mode should be also equal to zero 
in the frequency band studied, due to the fact that the cutoff frequency of the first higher 
order mode in the narrow section is / = cq/2w = 22.8 kHz. However we see a propagative 
band in the transmission term T( 1)0 ),(i,o) appearing at frequencies slightly above the cutoff 
frequency of the first antisymmetric mode (1,0), / = Cq/2W = 7622 Hz, generated by the 
TRs. Figure 2B zooms in the frequency band containing these resonances. We observe 5 
transmission peaks (as many as coupled cavities in the waveguide) between 7700 Hz and 
7800 Hz, corresponding to the transmission of mode (1, 0) with amplitude equal to 1. In 
this figure, we also show the terms T(i ;0 ),( 3 ,o), Tpqjxs.o) and Tqq^qq), which indicate the 
coupling of the transmitted mode (1,0) with incident modes (3,0), (5,0) and (7,0). All 
these modes, except mode (1,0), are evanescent at these frequencies. However, Figure 2C 
puts in evidence the ability of the device to convert evanescent waves into propagative 
ones. Remarkably, we observe conversion of modes (3,0), (5,0) and (7,0) with ampli¬ 
tudes bigger than 1, demonstrating the possibility to transfer efficiently subwavelength 
information through the device. 

To test these ideas experimentally, we measure the transmitted amplitudes of the 
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Figure 2: Transmission properties of the imaging device. (A) Modulus of the transmission 
coefficients T( 0 0 ),(o,o) (blue line) and T( 1)0 ) : ( lj0 ) (black line). The propagating band around 
7.75 kHz is due to the excitation of the trapped resonances in the waveguide. (B) Zoom 
in of the band containing the trapped resonances. (C) Measured transmitted amplitudes 
of modes (0,0) and (1,0). The blue solid line and the black dashed line represent, re¬ 
spectively, the amplitudes of modes (0,1) and (1,0), measured without object. The black 
solid line represents the amplitude of mode (0, 0) when an object is placed close to the 
waveguide input. 
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plane mode, ^4(o,o)> an d the first two antisymmetric modes ^4(o,i) and ^4(i,o) through the 
waveguide shown in Figure IB. The waveguide was placed vertically inside a box lined with 
absorbing foam, and a 22 mm diameter loudspeaker (Clarion SRE 212H) was placed on the 
waveguide axis 280 mm away from the waveguide input section. The transmitted pressure 
was measured with 4 1/4 inch (6.35 mm) microphones (G.R.A.S. 40BD), placed flush with 
the inner wall of the waveguide, 20 mm away from the last narrow section, as indicated in 
Figure 1A. The waveguide output was Elled with absorbing foam to minimize backward 
reflections. The microphones were placed at the midpoint of each wall, coinciding with 
the nodal lines, in order to separate the contribution of each mode. The transmitted 
amplitudes are obtained as A( 0 ,o) = ( Pi + Ps)/2 = (p 2 +Pa)/2, A( 0 ,i) = (pi — pz)/2 and 
v4(i >0 ) = ( P 2 — P 4 )/ 2 , where pi to p 4 are the complex pressure measured by microphones 1 
to 4. The pressure was measured using phase-sensitive detection to minimize noise. 

The transmitted amplitudes of modes (0,0), blue solid line; and (1,0), black dashed 
line, are shown in Figure 2C. The plane mode amplitude exhibits a band gap in the band 
[3.7,11] kHz, as predicted in Figure 2A. Since the acoustic source is placed symmetrically 
with respect to the waveguide longitudinal axis, none of the antisymmetric modes are 
excited, and therefore the amplitude of mode (1, 0) is close to zero in the whole frequency 
range. This situation changes when an object is placed close to the waveguide input. 
The solid black line in Figure 2C represents the amplitude of mode (0,0), measured 
when the edge of an aluminium plate was placed in front of the waveguide input at 
approximately 1 mm distance. In this configuration, the rapid variations of the acoustic 
held around the edge couple with the antisymmetric modes [Eq. (|2])], which in turn excite 
the TRs, and transmit signals around the TRs frequencies (see the inset of Figure 2C). In 
the experiments, we distinguish only 3 transmission peaks (instead of 5, as predicted in 
Figure 2B), and attribute this to inherent losses in the waveguide, not taken into account 
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in our model. 


The above results suggest that, at the TRs frequencies, the waveguide can be used 
to image acoustic field components with large perpendicular wavenumber. Let p J (y,z) 
be the acoustic field incident on the input waveguide cross-section S. The coupling of p 1 
with the waveguide modes is given by the prejection 

A (m,n) = f J^^m^dydz. ( 2 ) 

From Eq. ([2]), it follows that slow spatial oscillations of the incident field will couple mainly 
with the plane mode, which cannot propagate in the waveguide. In contrast, fast spatial 
asymmetric variations of the incident field will couple with the high-order antisymmetric 
modes and transmitted through the waveguide. However, since for any antisymmetric 
mode Eq. ([2]) vanishes if p 1 is symmetric, symmetric excitations cannot be transmitted. 

To test the imaging properties of our device, we performed a series of experiments 
to image the edges of different objects. Figure 3A shows a ID scan of the edge of an 
aluminum plate, the edges of a 32 mm wide aluminum plate are shown in Figure 3B, 
and the edges of a 10 mm wide aluminum rod are shown in Figure 3C. The frequency 
chosen is / =7740 Hz (A = 44 mm), corresponding to the maximum of A( 10 ) bi Figure 2C. 
The experimental results (solid red lines) were compared to finite elements simulations 
performed with Cornsol Multiphysics (black dashed lines). The single edge (Figure 3A) 
generates a sharp peak in the transmitted intensity of mode (1,0). The resolution, defined 
as the full width at half maximum of the peak (FWHM), is 0.22A. The two edges of the 32 
mm plate are seen as two narrow peaks (Figure 3B). We notice that these two peaks are 
slightly sharper (0.19A) than that of Fig. 3A. The reason for this is that when the position 
of the object is symmetrical with respect to the waveguide axis (y — 0), the transmitted 
intensity drops to zero, which shrinks the peaks generated by the edges. A special situation 
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arises when imaging small objects, whose edges are separated by a distance close to the 
device resolution (« 0.2A = 9 mm). Remarkably, the device still generates two sharp 
peaks with FWHM = 0.13A (Figure 3C). However, the resulting image does not represent 
the actual object size, but an object slightly larger. This is because the peaks maxima 
are shifted due to the intensity drop at y = 0. 



Figure 3: ID images of (A) an edge of an aluminum plate, (B) the two edges of a 32 mm 
wide aluminum plate, and (C) the two edges of 10 mm wide aluminum rod. The insets 
show the imaged objects, in which the dashed red line represents the scanned region. 


We have tested experimentally device’s ability to image 2D objects. In particular, 
we have imaged a 10 cm diameter plexiglas disc and the ETH Zurich logo, made of a 
rigid thermoplastic (Figure 4). Figure 4A shows the total transmitted intensity / = /(o,i) 
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+ /(i,o), when imaging the plexiglas disc (only the upper half disc is represented). The 
intensity is maximum at the edges of the disc (represented by the dotted line). The 
image also shows other features with lower amplitude, that can be generated by unwanted 
reflections in the experimental setup. A clearer image is obtained by reducing the dynamic 
range to one half of the maximum intensity, Figure 4B. A semicircle is clearly observed 
in the image, with FWHM « 0.2A. As mentioned above, an interesting aspect of this 
technique is the possibility to visualize edges aligned along different directions. This is 
achieved by visualizing the intensities 7(i ;0 ) (for horizontal edges) or 7( 0; i) (for vertical 
edges) separately. In Figure 4C, we visualize only 7(i,o) and observe a maximum in the 
region where the edge is horizontal, which vanishes smoothly as the edge becomes vertical. 
In Figure 4D, we visualize only 7(i,o) and observe that the intensity is maximum at both 
sides and vanishes as the edges become horizontal. 

Figures 4E-4H show images of the ETH Logo. The letters are 15 mm width (A/3) and 
the separation between letters varies between 10 mm (A/4.4) and 15 mm. We note that 
this situation is considerably more challenging than the previous cases, since the object 
contains a much larger amount of subwavelength information. The full dynamic range 
image, Figure 4A, shows intensity maxima coinciding with the edges of the letters. By 
reducing the dynamic range to one half, Figure 4F, the resulting image reveals most of the 
features of the object, except for the lower step of letter ,? T ,? . The edges of letters ”E” and 
” H” appear clearly in the image, as well as the upper part of letter ” T”. Figures 4G and 4H 
show the intensities 7(i,o) and /(o,i), which allow visualizing horizontal and vertical edges, 
respectively. These figures demonstrate that directional edge detection is also possible in 
this more complicated case. 

We foresee the ability to scale the fabrication of these devices to sizes of interest for 
ultrasonic imaging, to improve current visualization technologies in medical and non- 
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destructive evaluation applications. Moreover, since trapped modes also exist in elec¬ 
tromagnetic waveguides [2Si ES ], our results may suggest the design of analogous edge 
detection devices for optical waves. 
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Figure 4: (A) to (D), images of a 10 cm diameter disc. (E) to (H), Images of the ETH 
Zurich logo. A picture of these objects is displayed above these figures. (A) and (E) shows 
the normalized total intensity, / = /(i,o)+^(o,i)- (B) and (F) also show the total intensity 
but limiting the dynamic range to [0.5,1]. (C) and (G) represent /(i,o), which enables to 
visualize only horizontal edges. (D) and (H) represent which enables to visualize 

only vertical edges. 
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Supplementary Material 


SI Trapped resonances in waveguides with irregular, 
symmetric cross-section 


Trapped resonances exist in a variety of waveguide configurations (see, e.g ., (27) and 
references therein). In this section, we briefly describe the physical mechanism leading to 
the formation of trapped resonances in waveguides with irregular, symmetric cross-section, 
as the one shown in Figure 1A. 

Changes in the cross-section generate a modal coupling between modes propagating 
on both sides of the discontinuity. Due to the symmetry of the cross-section in y and 
z, symmetric modes on one side cannot couple to antisymmetric modes on the other 
side (their inner product vanishes). On the other hand, the cutoff frequency of the first 
antisymmetric mode in the wide segments, / = Co/2W, is below the cutoff frequency of 
the same mode in the narrow segments, / = Cq/2w. If this mode is excited in the wide 
section at a frequency between / = cq/2W and / = cq/2/w, it cannot couple to any 
propagative mode in the narrow section, and therefore it remains trapped. This situation 
is illustrated in Figure SI. 

evanescent antisymmetric modes propagative antisymmetric modes 
in the narrow section j n the narrow section 



0 1 —-— J k 

trapped 

resonances 


Figure SI: The cutoff frequency of the first antisymmetric mode in the narrow section 
is / = c 0 /2w (or k = tt/w) and that in the wide section is / = c 0 /2W (or k = tt/W). 
Trapped resonances appear in the band k G [n/W, tt/w[. 
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The TRs can be computed using finite elements by calculating the eigenmodes of 
the waveguide with perfectly matched layers (PML) on both extremities. The role of 
the PML is to take into account acoustic radiation. In a waveguide terminated with 
narrow segments, trapped resonances are represented by solutions with real resonance 
frequency, since the waveguide cannot radiate energy towards the extremities. Instead, in 
a waveguide terminated with wide segments, as the one considered here, the waveguide 
can radiate energy through the first antisymmetric mode, so that the TRs have a small 
imaginary part. 



Figure S2: Position of the first five trapped resonances in the complex A;—plane. The 
lowest images display the fields corresponding to the first and fifth resonance. These 
resonances generate the transmission peaks observed in Figures 2B and 2C. 

Figure S2 shows the first five TRs in the complex k— plane. These resonances are very 
close to the cutoff frequency of the first antisymmetric mode, / = Cq/2W (or k = 7r/IP), 
and generate the five sharp peaks observed in the transmission coefficients of Figures 2B 
and 2C. Their resonance frequencies (real part) are 7733 Hz, 7754 Hz, 7783 Hz, 7812 Hz 
and 7834 Hz, very close to the peaks observed in Figures 2B and 2D. Their quality factor, 
Q = Re{/c}/Im{fc}, is of the order 10 5 . 
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S2 Computation of the transmission matrix, T 


We have calculated the transmission matrix of the waveguide shown in Figure 1A using 
the multimodal method. Assuming perfectly rigid boundaries, the pressure held in the 
waveguide is the solution to the following problem (time dependence e Jujt is omitted): 

f (V 2 + k 2 )p(x, y, z), V(x,j/)efi, /^ 

\ d n p(x,y,z) = 0, V(x,y)eidQ, 

where k = cu/co is the wavenumber, O is the air inside the waveguide, dtt represent the 
boundaries and d n denotes the normal derivative with respect to the boundaries. 

The pressure held is developed on the basis of eigenmodes of the corresponding trans¬ 
verse section as 


p(x, y, z) — Y 4>( m , n ){y,z), (S4) 

m,n= 1 

with A( m n ) and 5( m , n ) the amplitude of the forward and backward modes, respectively, 
/3(m,n) — {k 2 —a 2 m n )) 1//2 the longitudinal wavenumbers, where Q?( m , n ) = [(mn/w) 2 + (■ nn/w) 2 ] lA 
in the narrow section and cc( min ) = [( rmr/W ) 2 + (mr /W) 2 ^^ 2 in the wide section are the 
traverse wavenumbers. The eigenfunctions <j>(m,n)(y, z ) are 

(j)(m,n)(y,z) = ~\/(2 — 8 m0 )(2 - 5 n0 )cos (y - ^)) cos (z - ^)) (S5) 

in the narrow section and 


0(m,n) (Zb %) 


1 

w 


V(2 — ^mo)(2 — S n o) 


COS 



in the wide section, where 5 is the Kroneker symbol. 


The scattering matrix of each element forming the waveguide (straight segments and 


discontinuities) is given by 

R t ,/ 

_ i i 

0i ~ [ T, R' 


(ST) 
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where i = 1,2,... ,1, with I = 13 the number of scattering elements, R, and T, are the 
reflection and transmission matrices for right-going incident waves, and R' and T' are 
the reflection and transmission matrices for left-going incident waves. 

The scattering matrices of the straight segments (i = odd) are given by 


S* 


[0] E 
E [0] 


(SB) 


where E is a diagonal matrix containing the terms e j/3( - m ’ n ) Ls , with L s the length of the 
segment, and [0] is the zero matrix. The scattering matrices of the discontinuities {i = 
even) are calculated from the continuity equations of pressure and normal velocity, given 
by 


pi 0 — pi 7 '') 


(S9) 


where superscripts (l) and (r) indicate quantities on the left and on the right of the 


discontinuity, respectively. Inserting Equation (S4) into Equations (S9) it is possible to 


obtain the scattering matrix for a sudden expansion as 

S,: = 


Ra T b 
T„ R h 


(S10) 


and the scattering matrix for a sudden narrowing as 

S ( = 


Ri, T a 
Tfe R a 


(Sll) 


with 


R 


I + F (Y (r) ) 1 (*FY W ) I * If (Y (r) ) 1 (*FY (/) ) - I 
T„=(yM) _ 1, FYW (I-Ra), 

R b = (Y (r) + ^Y^F)" 1 (Y (r ) - *FY (/) F) , 

T 6 = F(I + R 6 ), 


-1 r 


-1 
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where I is the identity matrix, F is the matching matrix, containing the inner product 
between the modal basis, {4^(mny$ik)l = I s ( $mn)$jk) < ^y < ^ z, > Y is a diagonal matrix 
containing the admittance terms and *F is the transpose of F. 

Defining the matrix operator * allowing assembling two scattering matrices, S* and 
Sj +1 , as (see Ref. (30)) 


Si*S 


i+1 — 


T i+1 (I - R(Ri+i) _ 1 T { 

R) + T'(I — R ?:+ 1 R') 1 R, i R' 


R' 


i+1 


+ Tj + i(l — R'Rj+i) 1 R(T ' +1 
T'(I - Ri + 1 Ri) _ 1 Ti +1 


(S12) 

The global scattering matrix of the waveguide, S, is calculated by assembling consecutively 
the scattering matrices of each element as 


S 


R T 
T R 


Si * S2 * • • • * S/, 


from which we extract the transmission matrix, T. 


(S13) 
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